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ABSTRACT 



Neglecting the ageostrophic circulation is not very 
realistic for ocean circulation in coastal waters and 
particularly in regions with strong temperature and salinity 
gradients like the Farallones National Marine Sanctuary. The 
C-Vector method has been developed to diagnose the three 
dimensional flow including ageostrophic circulation. This 
thesis uses the C-Vector method and data collected from the 
R/V Point Sur in the area of Farallones National Marine 
Sanctuary in order to describe and illustrate the three 
dimensional pseudo-vorticity fields of the ageostrophic 
circulation and the total circulation (geostophic and 
ageostrophic) of the region. 
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I . INTRODUCTION 



A. GULF OF THE FARALLONES NATIONAL MARINE SANCTUARY 
1. Geophysical Characteristics 

The boundaries of the Farallones shelf and slope 
region extend from Point Ano Nuevo on the south to Point Reyes 
on the north and westward from 10 km to 90 km off shore. 

This continental margin, which includes the coastal 
zone off San Francisco and is characterized by a shelf with a 
maximum width of 45 km, is the widest shelf along the entire 




Santa Cruz mountains (Shepard, 1941) . 

west coast of the United States (Figure 1) . This is, however, 
narrower than most of the shelves off the eastern and southern 
coasts of the United States. The shelf width reduces 
gradually to the north. Another characteristic that 
distinguishes the region is that the submarine canyons are 
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almost completely limited to the continental slope, beginning 
well out from the coast. 

West of Point Ano Nuevo, the direction of the 
continental slope veers northwest. The steepest part of the 
slope branches away from the shelf border north of latitude 
37° and a relatively smooth gradient extends to depths of 
1,000 to or 1,300 m. West of the Farallones Islands, the 
steepest zone is close to the shelf margin, but farther north 
and beyond Cordell Bank, there is a circular shape to the 
upper limit of the continental slope. 

The fact that most canyons are restricted to the 
continental slope in this region may have some bearing on the 
relatively smooth character of the upper parts of this slope 
(Figure 2) . However, there are many canyons on the outer 
slope. One of these canyons is the Pioneer Canyon between 
stations 8, 9 and 10 of section A and stations 11, 12 and 13 
of Section B. The maximum depth of this canyon is about 3,400 
m between stations 10 and 11 (Figure 3) . 

According to earlier studies (Stetson, 1936 and 
Shepard et al . 1941), currents in canyons might alternate in 
direction, but more often the strongest currents are down the 
canyon. It was thought in these early studies that the period 
of alternation was not tidal. 

Shepard et al . , (1979) believed that in areas with 
small tidal ranges, where short period alternations occur at 
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considerable depth, the tidal current oscillations are 
subordinate to regional current oscillations at other 
frequencies . 




Francisco (Shepard, 1941). 



An oval seamount, the Pioneer Seamount lies northwest 
of Pioneer Canyon and southwest of Farallones Islands. South 
of Pioneer Canyon is another oval seamount, Guide Seamount. 
Both rise about 1,000 m above the sea floor at the base of the 
slope, and each has more than one summit. These features may 
also play an important role in the water circulation and the 
currents of the region. 

2. Water Masses 

To study the water masses of the Farallones National 
Marine Sanctuary, the characteristic water masses in the North 
Pacific which may influence the currents west of San Francisco 
must be analyzed. 
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slope study. 

According to Sverdrup et al . (1964) the water mass 
which dominates the upper waters of the North Pacific is the 
subarctic water mass. The subarctic water mass is 
characterized as the water mass northward of latitude 45° 
North. For the sake of convenience, this name is used for the 
entire water mass above the latitude of 23° North. This water 
mass is advected eastward and then southward. 
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Between 22° and 45° and depths below 300 m, 
temperature and salinity (T-S) curves show waters intermediate 
in character between subarctic North Pacific and the 
Equatorial water. The water at the region of interest is 
formed by the mixing of these two large and well defined water 
masses. It was determined by Collins (1992) that a poleward 
penetration of Equatorial water (associated with the under 
currents) occurred along shore over the continental slope at 
the Gulf of Farallones Islands between 16 and 21 May 1991. 
Also, the waters were found to be warm and fresh showing 
characteristics of subarctic waters (at station 48) . 

Sverdrup and Fleming (1941) and Tibby (1941) 
discovered that near the coast (30 km) , the water mass had 
characteristics closer to the equatorial water mass than to 
subarctic water mass. It is obvious that close to the coast, 
the subarctic North Pacific water mass has less effect on than 
Equatorial Water mass characteristics. 

According to Hickey (1979) the California Undercurrent 
is a northward subsurface flow off the west coast of North 
America over the continental slope. The main flow, which 
carries equatorial type water, was found at 200 m depth. 
During the Farallones shelf and slope study, the undercurrent 
had maximum observed velocities at a depth of 100 to 200 m. 
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It has a jet-like structure, both vertically and horizontally, 
and extends to the bottom over the slope. 

The most permanent current of the shelf circulation 
during upwelling is the coastal jet which is a southward flow 
with a maximum speed located between 5 km and 25 km from the 
coast. Huyer et al . (1978) and Huyer et al . (1975) found that 
the seasonally averaged jet is stronger during the spring than 
in the summer. 

According to Yoshida (1967) , "processes near coastal 
boundaries respond greater and faster to time variable winds 
than those in the central oceans." During the first two weeks 
of September 1979, there were weak westerly winds which were 
not favorable for the maintenance of strong upwelling. At that 
period Breaker and Gillilard (1981) used IR images from the 
Advanced Very High Resolution Radiometer (AVHRR) to show the 
upwelled waters, near San Francisco, enter and leave San 
Francisco Bay. From the observed time sequences the upwelled 
waters appear to reflect the strong tidal influence at this 
location (Conomos, 1979) . 

The Davidson current which is a northward near-shore 
flow is commonly observed north of Point Conception during the 
winter months (Reid et al . , 1958), but can occur in other 
seasons under favorable wind conditions (Wyatt et al . , 1972) . 
The fact that drift bottles released north of Point Conception 
have been recovered as far north as Vancouver Island 
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(Schwartzlose, 1963; Crowe and Schwartzlose , 1972) is strong 
evidence that the flow is usually continuous alongshore. 

The California Current (CC) which is more fully 
developed during spring and early summer has its core 
approximately 100 - 200 km offshore (Chelton, 1984) and for 
this reason, does not greatly affect the area of study which 
extends only 90 km off shore. 

3 . Weather Conditions 

Weather conditions during the cruise are shown in 
Figure 4. The hourly averaged wind vectors are from west- 
northwest as usually happens during this season in this area 
(Huyer, 1983). The wind speed ranges from 3.8 to 17.2 m/sec 
and contributes to the upwelling process in the region. Air 
temperature is relatively low, between 9.8° C and 11.6° C. 

B. GENERAL CLIMATOLOGICAL PICTURE 

The data used in this thesis was gathered as part of an 
Environmental Protection Agency (EPA) funded plan to study 
oceanographic currents and hydrographic conditions in the area 
of the Gulf of the Farallones National Marine Sanctuary 
(Figure 3 ) . 

According to the schedule, five surveys were conducted. 
This survey was carried out between May 16 to 21, 1991 aboard 
the research vessel Point Sur and was the second survey of the 
five . 
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board R/V Point Sur during the cruise of May 16 - 21, 1991 

(Jessen et al . , 1992). 

During this cruise, a three-dimensional map of 
hydrographic data was obtained in order to better understand 
the currents of the area. 

Additionally, data was acquired from Advanced Very High 
Resolution Radiometer (AVHRR or AVHRR/2) instruments on board 
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TIROS -N/NOAA polar orbiting satellite and interpreted as Sea 
Surface Temperature (SST) images. 

Table I. ISLANDS, ROCKS AND SHOALS IN THE GULF OF THE 
FARALLONES NATIONAL MARINE SANCTUARY. 



Islands , Rocks 


Latitude 


Longitude 




and Shoals 




North 




West 




1 . S . E . Farallon Is. 


37 


42 


123 


00 




2. Mid. Farallon Is. 


37 


44 


123 


0# 




3. North Farallon Is. 


37 


46 


123 


0® 




4 . Fanny Shoal 


37 


47.5 


123 


10 . 5 




5. Noonday Rock 


37 


47.5 


123 


• 0 




6. Hurst shoal 


37 


41.5 


123 


00 


1 

j 


7. Maintop Is. 


37 


41.8 


123 


00 





The locations of the Farallones Islands, rocks and shoals 
near stations 22, 23 of Section C and 37, 38 of Section D, are 
listed in Table I. 
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C. PURPOSE OP THE STUDY - CONTENTS 



This research investigates the three-dimensional 
ageostrophic pseudo- vorticity fields and the pseudo- vorticity 
fields of the total circulation near continental shelves and 
slopes extracting the necessary information from the 
hydrographic and wind data sets (Figure 5) . 

The California Current (CC) , which is the eastern part of 
the subtropical North Pacific gyre and extends from Washington 
State to Baja California, is the main current at the surface. 
The California Undercurrent or Inshore Current (IC) with 
poleward flow is the main current below the surface in this 
region. The geostrophic part of these currents can be 
directly obtained from the hydrographic data sets (CTD data 
sets) . 

The ageostrophic portion of these currents can be 
indirectly obtained by the C-Vector which can be calculated 
from hydrographic data and surface wind stress data. 

In this thesis, Chapter II covers the C-Vector concept and 
the techniques for the C-Vector calculations. Chapter III 
describes instruments, methods and data processing techniques. 
Chapter IV provides detailed analysis of the plots of the C- 
Vector for the ageostrophic circulation and the C-Vector for 
the total circulation (geostrophic and ageostrophic) . 
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Figure 5. Schematic diagram for the calculation of 
ageostrophic pseudo-vorticity fields and then the pseudo- 
vorticity fields for the total circulation. 



Finally, Chapter V summarizes the results and offers 
recommendations for future work. 
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II. C- VECTOR CONCEPT AND TECHNIQUES 



A. OBJECTIVE 

The main objective of the C-Vector concept is to compute 
a three-dimensional pseudo- vorticity fields from the density 
fields and its derivatives, and then to obtain a three- 
dimensional circulation. 

B . THEORY 

1. Meteorology 

The general idea was first proposed for atmospheric 
circulation with the Q-vector concept of Hoskins et al . 
(1978). Hoskins et al . found that the Q-Vector was obtained 
from the w-equation and can be explained on the synoptician's 
chart more precisely than in earlier studies (Sutcliffe, 
1947) . This theory is very helpful in diagnosing the synoptic 
and frontal vertical circulation. A disadvantage of this 
concept is that it ignores the barotropic part of the 
rotational ageostrophic flow. According to Xu (1992) this 
problem can be overcome by the C-Vector concept which resolves 
the ageostrophic circulation more conveniently and allows a 
new way of solving the C-Vector equation. The C-Vector 
equation is a three-dimensional extension of the Q-Vector 
equation (QG diagnostic equation) . 
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Figure 6. Vertical Temoerature and Salinity plots (Jessen at 
al. 1992) . 
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2. Physical Oceanography 

The C-Vector concept can also be applied to oceanic 
circulation (Chu, 1992) . With the help of this concept and a 
set of hydrographic and surface wind data, the three- 
dimensional pseudo-vorticity fields of ageostrophic 
circulation in the oceans are diagnosed. Also revealed are the 
three dimensional pseudo-vorticity fields of the total 
circulation which includes both geostrophic and ageostrophic 
flow. 

Only CTD data are needed for this calculation of the 
currents, obviating a need for more expensive and difficult 
direct measurements of ocean currents in sufficient quantity. 

There is never a perfect geostrophic balance. A 
geostrophic balance is traditionally assumed for calculation 
of the currents because the ageostrophic part is considered to 
be very small. However this assumption is of questionable 
validity especially for coastal waters of the Farallones 
National Marine Sanctuary where there are significant 
temperature and salinity gradients like the one at Section E 
and D (Figure 8 ) . 

The thermal wind equations, 



— (p fv) = -g 
dz M 9 y dx 



-$iW% 



(1) 
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(where p is the sea water density, f is the Coriolis parameter 
jf=2Qsin0) are used for the calculation of the geostrophic 
currents (U g ,V g ). The geostrophic currents are not, however, 

f 

the total currents. In order to have more realistic flow 
fields the ageostrophic circulation must be included. 
Following Chu (1992), the basic equations for coastal waters 
including ageostrophic circulation, wind stress and excluding 
buoyancy are: 






( 2 ) 



<^.V)v + f(u-^)=if 



(3) 



( 



-JL + VV) b+N z w=0 
dt 



(4) 



du + dv + dw 
dx dy dz 



(5) 



where N is the Brunt -Vaisala frequency (assuming N = 
constant) , f is the Coriolis parameter (on an f - plane) , 
V=(u,v,w), and (Y y ,Y x ) is the turbulent momentum flux due to 
surface wind stress (r x ,T y ), 



pc | 





(6) 
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The decomposition of the total current velocity into 
geostrophic velocity and ageostrophic circulation results in: 

V=V v *V a (7) 



Using the above equation (7) and after the 
decomposition of the basic equations (2) to (4) , the system of 
equations will have the form, 

t K(V a ) =3t(V g ,p,Yy,Y x ) (8) 



where %, 9? are specific differential operators. From 
functions !K and 9?, the three-dimensional ageostrophic 
circulation can be calculated using equation (8) together with 
CTD and surface wind data. 

3. Development of C-Vector equations (a Diagnostic 
System) 

Equations (2) to (4) , for coastal waters on an f - 
plane, can be decomposed using equation (7) . The result will 
be the following equations: 



-fv =$Xl-(JL + v -V) U 
a dz dt 9 ’ 9 



(9) 



fu 



a 



dY y 

dz 



( — + V -V) V 
dt 9 9 



(10) 
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( 11 ) 



N 2 W a = -(-^+V g V)b 



du„ dv , dw . 

+ ^-£ +^-£ =0 



( 12 ) 



dx dy dz 

Taking Che derivative of equation (9) with respect to 
z, and f times the derivative of equation (11) with respect to 
y, and adding them together (similarly working for equations 
(10) and (11)) results in Che ageostrophic pseudo- vorticity 
equations : 



~ + 4~ ( N 2 W ) =2C 1 

dz a dy a 



( 13 ) 






( 14 ) 



4 - - 4 - ( f2u a) =2C 3 

dx a dy a 3 



( 15 ) 



where : 



du„dv, a v-.acr , t yr- 

1 dy dz dy dz 2 dz 2 



( 16 ) 



Q — f(- 



dU Q dV g _ dV„ dU Q , , f d z Y y 



) + 4 - 



dz dx dz dx 2 dz 2 



( 17 ) 



3 dx dy dx dy 2 dz dx dy 



( 18 ) 
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C a s(C lf C 2 f C 3 ) 



( 19 ) 



The above equations yield the three components of the 



C a - Vector. Having the C a - Vector, the three dimensional 
ageostrophic circulation can be found. 



for the total flow, geostrophic and ageostrophic. The three 
components of this vector can be found with the following 
process . 

Using equation (13) for the total flow (v=v g +v a and 
w=w g +w a but w g =0 => w=w a ) : 



Another important vector is the C-Vector (C x ,C y ,C z ) 




*~4z < f2l V +2C 1 *2C x f 2 - 



( 21 ) 




C i _ 1 dv <? 



( 22 ) 



f 2 2 dz 



Working similarly, the other components can be found: 



Cy f2 2 dz 




( 23 ) 



Cz f 2 2 dx dy 




( 24 ) 
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4. Physical Significance of C-Vector 

There are two important physical meanings of C-Vector. 
a. Direction of C-Vector 

There is a non-dimensional formula which can 
appropriately show the relation between C-Vector and 
ageostrophic pseudo-vorticity : 



VxV a =2R 0 C (25) 

V-V a =0 (26) 

where R 0 is the Rossby number (R 0 =U/ ( fL) ) • 

Equation (25) shows that C-Vector is proportional 
to the ageostrophic vorticity. A C-Vector stream line can be 
considered as an ageostrophic vortex line (Figure 7) . 

b. Relation between w a and vertical vorticity of C- 
Vec tor. 

This relation can be found from: 
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Figure 7. Ageostrophic circulation induced by C-Vector in 
three-dimensional space. 

-ic-Vx ( equation 26) — 

V 2 iv a = -2i? 0 ic , VxC (28) 

Equation (28) demonstrates that the vertical 
velocity (w a ) is caused by the vertical vorticity of C-Vector 
(Figure 8 ) . 

Define : 
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Figure 8. Vertical velocity w a induced by horizontal 
rotation (vorticity) of C-Vector. 



W=2(^-^) ( 29 ) 

ox oy 

Combining equation (29) with equation (28) results 
in the following: 



V^w a =-R 0 W 



( 30 ) 



Equation (30) is the Poisson equation. If normal 
mode solutions are considered, the vertical velocity w a has 
the same sign as ^ , therefore ^ will be very useful because 
its positive (negative) values show positive (negative) 
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vertical ageostrophic velocities. The solutions 
information regarding upward (upwelling) and 
(downwelling) motion. 



provide 

downward 
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III. DATA COLLECTION AND PROCESSING 



The R/V Point Sur departed from the Moss Landing Harbor on 
May 16, 1991 with a combined scientific crew from Naval 

Postgraduate School (NPS) , Moss Landing Marine Laboratory 
(MLML) , Scientific Applications International Corporation 
(SAIC) , California Academy of Science (CAS) and the 
Environmental Protection Agency (EPA) . The purpose of the 
cruise was to study oceanographic currents and hydrographic 
conditions in the area of the Gulf of the Farallones National 
Marine Sanctuary. The ship proceeded as planned along five 
across-shore transactions 20 km apart beginning with CTD 
station 1 (Figure 3) . A total of 50 CTD casts were conducted 
to within approximately 25 m of the bottom. A brief ADCP 
survey was conducted as well (Figure 6) . 

This data was used to determine the pseudo- vorticity 
circulation of the region using the C-Vector concept. 

A. DATA COLLECTION 
1. CTD 

All CTD stations with their positions, weather 
conditions and maximum depth (m) are shown in Table II. 
Stations 1, 20, 21, and 40 were not included because they were 
located outside the grid (Figure 9) which was used for the 
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Table II.CTD STATION LOCATION, WEATHER AND MAX. DEPTH 



Stat. 


Lat. 


Lat. 


Lon. 


Lon. 


Wind 


Wind 


Air 


Max. 


No. 


Deg. 


Min. 


Deg. 


Min. 


Dir. 


(m/s) 


Temp. 


Depth 




N 


N 


w 


w 






(C) 


(m) 


19 


37 


5 9 


123 


18.1 


288 


7 1 


116 


2858 


9 


37 


10.2 


123 


S.S 


28f 


11.7 


11.3 


1440 


8 


37 


14 2 


123 


0.9 


291 


12 4 


11.3 


594 


7 


37 


16.3 


122 


56.4 


276 


11.5 


11.2 


390 


6 


37 


ie * 


122 


49 1 


289 


12.3 


10.7 


210 


5 


37 


10 s 


122 


49.1 


288 


8 8 


10.3 


110 


4 


37 


21 1 


122 


46.4 


356 


8.8 


10.2 


86 


3 


37 


23 8 


122 


40.7 


289 


7 7 


10.2 


86 


-> 


37 


25.1 


122 


36.2 


285 


10 8 


10.3 


52 


19 


37 


36 2 


122 


42.0 


301 


17.2 


10.0 


42 


18 


37 


3S 4 


122 


47.4 


307 


15 3 


10 0 


66 


19 


37 


33.6 


122 


53.0 


297 


14.6 


9.8 


92 


16 


37 


29 2 


122 


55.9 


289 


16.7 


9 8 


126 


15 


37 


27.9 


122 


59.1 


306 


15.4 


9.9 


342 


14 


37 


25.5 


123 


3 4 


295 


15.5 


10 0 


594 


If 


37 


23.7 


123 


7 . 1 


303 


13.9 


10 0 


806 


12 


37 


19 3 


123 


10.8 


303 


13.9 


10 • 


1400 


11 


37 


15.5 


123 


25.0 


316 


ii.T 


10 0 


2738 


30 


37 


24.6 


123 


32.5 


301 


11.3 


10.0 


2840 


49 


37 


29.0 


123 


23.4 


305 


12.0 


10.3 




28 


37 


33.2 


123 


14.8 


303 


10.8 


9.1 


1630 


27 


37 


35.4 


123 


10.1 


297 


8 S 


9.1 


1214 


26 


37 


37.3 


123 


6.4 


313 


7.7 


9.G 


714 


21 


37 


38.6 


123 


3 4 


289 


4 4 


10 • 


112 


24 


37 


40 .0 


123 


0.5 


289 


4 8 


10.0 


S6 


23 


37 


42.0 


122 


54 8 


286 


5.7 


10.1 


36 




37 


45.5 


122 


49 l 


286 


6.7 


10.1 


46 


33 


37 


53 8 


122 


53.3 


289 


7.2 


10 0 


S6 


38 


37 


50.9 


123 


3.3 


289 


6.7 


10.0 


76 


47 


37 


48.0 


123 


8.3 


293 


72 


9.8 


86 


36 


37 


46.7 


123 


9.2 


293 


7.4 


9.8 


16 


3S 


37 


45.3 


123 


12 0 


294 


6.4 


9.8 


366 


34 


37 


43 3 


123 


10 3 


297 


6.4 


10.0 


1054 


33 


37 


41.0 


123 


20.8 


292 


8 8 


10 0 


1674 


32 


37 


38 4 


123 


26 6 


295 


7.7 


10.2 


2712 


31 


37 


34.9 


123 


34.9 


313 


y.4 


10.7 


3050 


49 


37 


41.5 


123 


40 8 


305 


21 


10 0 


3146 


48 


37 


45.6 


123 


32.2 


289 


7.0 


10.5 


2412 


47 


37 


47.5 


123 


27.6 


299 


8.8 


10.7 


1620 


46 


37 


49.5 


123 


23 3 


298 


7.3 


10.5 


582 


45 


37 


50.9 


123 


20.2 


299 


3.8 


10.6 


112 


44 


37 


52.3 


123 


17.5 


289 


4.0 


10.7 


98 


43 


37 


55.0 


123 


11.7 


288 


7.2 


10.7 


90 


42 


37 


57.8 


123 


6.0 


275 


6.2 


10.5 


74 


41 


37 


59.1 


123 


3.2 


270 


4.7 


10.5 


40 
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calculation of the C-Vector. The CTD casts were made with a 
Neil Brown Mark III-B CTD attached to a General Oceanic 
rosette sampler for in situ water sampling. This model CTD is 
considered to be accurate to within ± 0.005 mmho, ± 0.005° C, 
and ± 6.5 db. Conductivity, temperature and pressure 
resolutions are 0.001 mmho, 0.0005° C and 0 . 1 db respectively. 

The CTD sensors were calibrated prior to the cruise. 
Coefficients calculated for these calibrations were applied to 
the raw data. Salinities were further calibrated by adjusting 
the calculated salinity to the salinity of water samples 
collected during the cruise. Further details of instrument 
calibrations can be found in the NPS data report for this 
cruise (Jessen et al . , 1992). 

2 . ADCP 

Continuous relative current velocity data was gathered 
aboard the R/V Point Sur using RD Instruments Vessel Mounted 
Acoustic Doppler Current Profile (VM-ADCP) , equipped with a 
four beam JANUS array, operating at a frequency of 150 khz. 

3 . Wind data 

During the cruise, an R. M. Young anemometer on the 
R/V Point Sur obtained 30 second averages of wind speed and 
direction (Table II) . 

4. Satellite 

The AVHRR Instrument on NOAA-11 satellite measured 
emitted radiation of the study area in five wavelength bands 
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(three thermal Infra Red (IR) Bands, one Near Infra Red (NIR) 
and one Visible Band) . 

B. DATA PROCESSING 

1. CTD 

The raw CTD data was processed on a PC using an EG&G 
Marine Instrument software package called CTDPOST. This 
program flags bad data values using a first difference 
criterion and allows the user to examine and interpolate 
across bad data if necessary. After the correction of the 
data, a dual time lag filter was applied to remove time lag 
spikes, and the data was averaged to 2 db. 

2 . ADCP 

The conversion of recorded ADCP relative velocity data 
into absolute velocity data took place on board ship. This 
process depends significantly on the accuracy of the ship's 
navigation instruments. The ship's position data was obtained 
from GPS and was also recorded on computer disks for later 
processing. 

Variability of the ship's speed, direction, and data 
collection interval, play significant roles in the ADCP 
accuracy and reliability. More details concerning ADCP data 
processing techniques can be found in Kosro (1985) . According 
to Kosro, accuracies of 4 - 5 cm/sec in the athwartship 
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component and 2-4 cm/sec in the fore -aft component are to be 
expected . 

3 . Wind data 

The wind data, collected underway, was stored on the 
hard disk of an HP310 computer. This data was transferred to 
5.25 inch diskettes upon return to port and computer processed 
to obtain hourly averaged wind vectors (Figure 4) . 

4. Satellite data 

The Satellite data was first calibrated, transferred 
to earth coordinates and temperature calculated pixel by 
pixel. All images were coregistered to a similar coordinate 
map of the study area. This work took place at the Naval 
Postgraduate School IDEA Lab. 

5. Geostrophic velocities 

All densities and geostrophic velocities were obtained 
from the calibrated CTD data using computer programs developed 
in the NPS oceanography department based on the equation of 
state of seawater (Fofonoff, 1985) . 

The u and v components are the across shore and along 
shore components of the geostrophic velocities. These 
components at a particular station were calculated using 
stations on either side of it. For example, utilizing 
stations 26 and 24, the v- component of geostrophic velocity of 
station 25 can be determined, and from stations 16 and 36, the 
u- component of station 25 (Figure 9) . 
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This does not work for the external stations. In this 
case the u and v components of the external stations were 
calculated using the external stations and the stations close 
to them. For example, for the determination of the u- 
component of geostrophic velocity at station 3, stations 3 and 
13 were used. 
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6. Reference level or Level of No Motion (LNM) 

The calculation of geostrophic velocities depends on 
the right choice of a reference level or Level of No Motion 
(LNM) . 

The hydrographic and ADCP data from the NPS Farallones 
Shelf and Slope study on May 16 to May 21, 1992 were examined 
carefully. It was determined from the plots of the vertical 
sections of 5 km averaged, across and along- transect ADCP 
velocities (for sections A, B, C, D and E) that the existence 
of a common horizontal reference level or LNM was uncertain. 
It was obvious that the reference level could not be 
horizontal over a great distance. 

In this study a large number of stations were located 
in shallow waters so the area was divided into two parts 
(Figure 10) . A shallow part which includes stations with 
depths less than 1,000 meters (over the continental shelf), 
and the deep part which includes stations with depths more 
than 1,000 meters (over the continental slope) . It was 
assumed that the reference surface was horizontal within the 
deep waters. 

According to Soule (1939) the water density field has 
to be analyzed first and the horizontal reference surface has 
to be placed in the layer where density is least disturbed. 

In analyzing the distribution of density differences 
between neighboring stations the depths were found where 
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Figure 10. A 3-D view of the stations over the shallow and 
deep waters. 

density differences changes sign or became zero. These depths 
represent the levels of minimum current shear. 

Figure 11, which illustrates the plots of the density 
differences between neighboring stations in the shallow part 
of Section D, shows that there is no common depth where the 
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Figure 11. Density differences between pairs of stations of 
Section D. 



density differences are zero. The choice of the LNM in 
shallow waters is therefore uncertain, and the reference layer 
was estimated to be the lower level of measurements between 
neighboring stations. 

The analysis of density difference distribution for 
deep stations offers better results. Figure 12 includes all 
the pairs of deep stations. A common depth of "zero" density 
differences is found at 1,400 meters. This level (1,400 m) 
was used as a reference level for deep waters. 
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There is another method of determining the LNM similar 
to the above process (Mamaev, 1955). The difference is that 
the curves for the vertical distribution of differences of 
specific volume at two neighboring stations are analyzed, 
instead of the differences of density. 

7. Computation of C-vector 

a. Calculation of C-vector (C lt C 2 , C 3 ) from the 

geos trophic currents. 

The computation was done using equations (16) - 

( 18 ) , where : 

• U g and V g can be calculated for all stations. 

• Coriolis parameter was assumed constant, f=2Qcos0 (where 
Q = 7.29 10' 5 sec' 1 and an average latitude 4> = 37° 30' 
North is used) 

• Y y and Y x , the turbulent momentum fluxes due to surface 
wind stress (r x ,T y ), are given by equations (6) 

Taking the average density of all stations it was 
determined that an average ocean density of the region was p 0 = 
1.026 gr/cm 3 . 

Wind stress (r x ,r y ) is obtained from : 



. t y =p«c 0 KK (31) 

where p a = 1.293 10' 3 gr/cm 3 (air density), C D = 1.3 10' 3 

(drag coefficient) , |V W | is the hourly averaged wind speed at 



33 



height 10 m and u w , v w are the x and y - components of wind 
speed respectively. 

The mixed layer depth (MLD) of the study area is 
needed to determine to what depth the wind mixing influences 
the ocean stratification. Observing all MLD from temperature 
plots during this study, it was determined that an average MLD 
(of all stations) is 15 m. This means that wind stress mixes 
the waters down to 15 m. 



A linear reduction of wind stress with respect to z 
was assumed (Figure 13) for the purposes of C-Vector 
calculation. 




depth . 

Jb. Computation of C-vector (C x , C y , C z ) due to the 
total flow (geos trophic and ageos trophic) . 

This computation depends on the equations (22) - 

(24) where C 1; C 2 and C 3 and geostrophic velocities are 
calculated as described above. 
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Having C x , and C y it is then easy to compute ¥ from 
equation (29) . Contours of C-Vector for ageostrophic flow, C- 
Vector for the total flow and are analyzed and compared in 
the next chapter. These three parameters are the main tools 
for the identification of the ageostrophic and the total 
(geostrophic and ageostrophic) pseudo-vorticity fields. 



35 



IV ANALYSIS OF PSEUDO -VORTICITY FIELDS 



A. DESCRIPTION OF THE ANALYSIS 

This chapter presents a detailed description of the 
pseudo- vorticity fields for the ageostrophic circulation 
(C 1 ,C 2 ,C 3 ) in the Gulf of the Farallones. The results of this 
description enable us to proceed further in the analysis of 
the pseudo-vorticity fields for the total geostrophic and 
ageostrophic circulation (C x ,C y ,C z ) . Then the vorticity of C- 
Vector (Sf r ) will be used to indicate the upward and downward 
motions in the area, and it will be compared to a satellite 
picture of the same area and the vertical ADCP velocities. 

1. Pseudo-vorticity fields ( c i' C 2 ,C 3) for the 
ageostrophic circulation 

a. Analysis of the x-component of ageostrophic 
pseudo-vorticity field 

C x is the ageostrophic pseudo-vorticity component 
of C-Vector in the x-direction or cross-shelf direction. This 
direction approximately indicates the cross polar vertical 
circulation pattern. The values for C 1 /f 3 have been plotted 
on the y- direction or along-shore direction. 

From Figure 9, notice that there are nine columns 
on the y- direction. Columns 1, 2, 3, 4 and 5 correspond to 
shallow water and are going to be analyzed first. 
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Figures 14, 15 and 16 show the plots of the first 
five columns. The arrows on the plots indicate the direction 
of rotation of the ageostrophic pseudo-vorticities and not the 
directions of currents. The cells with positive (negative) 
values of the ageostrophic pseudo-vorticity are turning 
counter-clockwise (clockwise) or cyclonic (anti-cyclonic). The 
direction along the positive (negative) axis is upward 
(downward) out of the page (into the page) . 

From these five plots it can seen that the 
ageostrophic pseudo-vorticities are increasing going offshore. 
At column 1, the maximum values of C 1 /f r3 are about +/- 0.1 and 
at column 5 the values have increased to 1.3. The first two 
columns have the negative ageostrophic pseudo-vorticity cells 
north and the positive cells south. The opposite occurs at 
columns 3 and 4 and at column 5 the positive ageostrophic 
pseudo-vorticity predominates from South to North. 

At column 5 the strong south positive ageostrophic 
pseudo-vorticity cell coincides with the area of large 
temperature and salinity gradients of Section A (Figure 6) . 
Also, the strong middle positive ageostrophic pseudo-vorticity 
cell coincides with the significant salinity gradient of 
Section D (Figure 6) . Column 5 has the largest positive value 
of the ageostrophic pseudo-vorticity fields. 

Figures 16, 17 and 18 show the four plots of 
columns 6, 7, 8 and 9 which lay over the deep waters. Column 
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7 shows a strong positive ageostrophic pseudo- vorticity field 
that dominates the area. Columns 6 and 8 have their negative 
ageostrophic pseudo-vorticity cells on the north side and 
Column 6 has one positive cell south. Column 8, however, has 
a positive cell followed by another negative cell south of the 
north positive cell. Column 9 has a significant negative 
north ageostrophic pseudo-vorticity field. Its position 
coincides with the large temperature and salinity gradient of 
the west side of Section E (Figure 6) . 

The significant temperature and salinity gradients 
are probably due to subarctic water underlying a surface lens 
of fresh and warm water (Collins, 1992). Figure 19 presents 
a better picture of the ageostrophic pseudo-vorticity fields 
through a 3-D plot of columns 9,7 and 5. The red lines are the 
solid positive values and the blue lines are the dashed 
negative values. 

The majority of the ageostrophic pseudo-vorticity 
fields are located at the upper 200 m. Waters below 200 m 
were found to have very weak ageostrophic circulation, so the 
significant ageostrophic circulation is close to the surface. 

b. Analysis of the y-component of the ageostrophic 

pseudo-vorticity fields 

The values for C 2 / f 3 have been plotted at the five 
Sections A, B, C, D, and E (Figures 20, 21 and 22) . Section 
A shows a relatively weak ageostrophic pseudo-vorticity field. 
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The plots of Sections B and C include four cells. Moving off 
shore alternating values of negative and positive cells are 
observed. Sections D and E have the highest values of the 
non-dimensional component C 2 /f 3 . The west (left) side of 
Section E has the maximum negative value of the ageostrophic 
pseudo- vorticity fields. This location coincides with the 
large off - shore temperature and salinity gradient of section 
E (Figure 6) . The values of C 2 / f 3 below 200 m are very small 
as can be seen from the plot of Section C (Figure 21) . 

The magnitudes of C 2 /f 3 (maximum value -7) are 
larger than the magnitudes C x / f 3 (maximum value 1.3). This 
probably happens because the average distance between the five 
sections (dy = 20 km) is bigger than the average distance 
between the nine columns (dx = 8 km) , and the derivatives with 
respect to y are smaller than the derivatives with respect to 
x. 

Figure 23 offers a three dimensional view of the 
ageostrophic pseudo-vorticity cells of Sections A, C and E. 
All the sections show that the ageostrophic circulation is 
located close to the surface and close to the locations of the 
large temperature and salinity gradients. 

c. Analysis of the z- component of ageostrophic 

pseudo-vorticity field 

As expected the values of vertical component C 3 /f 3 
are smaller than the values of the horizontal components 
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because the vertical ageostrophic vorticity is much smaller 
than the vertical geostrophic vorticity. 

The upper plot of Figure 24 shows C 2 /f 3 at the 
surface, without the wind affect. Notice that there are three 
negative ageostrophic pseudo- vorticity cells and one positive. 
The large values of C 2 /f 3 are concentrated at the north-west 
side of the area where the large temperature and salinity 
gradients were observed before. The lower plot of Figure 24 
shows the contours of C 2 /f 3 including the wind. The majority 
of the values of C 2 / f 3 are negative at surface. This probably 
happens because the area with the negative values coincide 
with the area of strong winds. A comparison of the plot of 
C 3 /f 3 at the surface with wind and without wind shows that 
without wind there is one cell with relative large positive 
values of C 2 / f 3 on the North-West side. 

The plot which includes the wind shows that the 
negative values predominate in the region, and the positive 
cell that existed without the wind has become smaller and 
weaker (the magnitude is smaller) . For C 3 / f 3 at 50 m (Figure 
25) , the three negative cells and the one positive cell exist 
at the same positions but with smaller values. 

At 200 m (Figure 25) , the three negative cells 
still exist almost at the same location but with smaller 
magnitudes. The previously formed positive cell is very weak, 
and a new positive cell appears at the south-east corner of 
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the area. Figure 26 shows a three-dimensional picture of the 
three plots of C 3 /f 3 at the surface and the depths of 50 m and 
200 m. 

Figure 27 shows a 3-D picture of C-Vectors of 
columns 9, 7 and 5. The large negative values of C- Vector 

(for ageostrophic circulation) were found at the west side of 
Section E where the large temperature and salinity gradients 
were observed. 

2. Pseudo- vorticity fields for the total circulation 
(C x ,C y ,C z ) ageostrophic and geostrophic 
a. Analysis of the x- component of C- Vector (C x ) of the 

pseudo -vorticity field 

The columns 1, 2, 3, 4 and 5 are located over the 
shallow waters. At these columns the negative pseudo- 
vorticity fields are dominant (Figures 28, 29 and 30) . The 
strongest cells are located south and north of column 5, and 
the south one coincides with the large temperature and 
salinity gradients at the middle of Section A. 

At the deep waters, the negative cells of pseudo- 
vorticity are dominant also (Figures 30, 31 and 32). At the 
north part of column 8 and 9, the maximum values of C x / f are 
observed (C x /f = -35) . In the same area, maximum values of 
C x /f 3 and C 2 /f 3 and C 3 /ir 3 were seen before. This is the area 
with maximum temperature and salinity gradients. The arrows 
show an upward flow at column 9. At column 8, there is an 
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upward motion south and north and at the middle downward flow. 
The arrows indicate the upward and downward motion of all 
columns . 

Figure 34 is a three-dimensional plot of C x /f for 
columns 9, 7 and 5 where the upward- downwards motion can be 
seen more clearly in three-dimensional space. 

b. Analysis of the y- component of C -Vector (C y ) of th e 

pseudo -vorticity field 

Positive pseudo- vorticity cells of C y dominate the 
five sections (Figures 34, 35 and 36) . Notice that a negative 
pseudo-vorticity cell exists at the west side of all Sections. 
The arrows show an upward motion off shore of all Sections 
followed by a downward motion (going east) . At Sections A, B, 
and C, there is an indication of upward motion at the east 
side of these sections. Sections D and E have the strongest 
pseudo-vorticity cells at their west side (maximum C y / f = 13 
at Section D) . This area contains the large horizontal 
temperature and salinity gradients. Figure 37 contains a 
better picture (3-D) of the pseudo-vorticity fields of 
Sections A,C and E. 

c. Plots of ¥ function 

Figure 3 8 shows ^ at the surface. The large 
positive values of ^ mark an upward motion at the west side 
of the region. In the same figure, there is an outer area 
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west of the Farallones Islands where there are positive values 
of ^ and these are probably upwelling waters. 

Between these two areas, there is an area with 
negative values of 'P (downward motion) . This area lies over 
the continental slope. A weak upward motion (positive values 
of 'i') appears at the east side of the region and south-east 
of Farallones Islands. At 20 m, the picture is similar to the 
surface . 

At 50 m (Figure 39) , the locations of upward and 
downward regions are almost the same as those seen at the 
surface and 20 m. At 200 m, the situation starts to invert, 
and the regions of upward motions become regions of downward 
motions. This means that the upward -downward motions do not 
extend as deep as 200 meters. 

Figure 40 is a 3-D plot of ^ at surface, 50 meters 
and 200 meters. In this color picture the solid positive 
values are red (upward motion) and the dash negative values 
are blue (downward motion) . 

B. VERIFICATION OF C-VECTOR VORTEX LINES 

1. Direct verification of C-Vector vortex lines by 'P 
function 

An important verification of the upward and downward 
water motions at Gulf of Farallones Islands induced by C- 
Vector can be achieved by the ^ function. The positive values 
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of ¥ give the upward motion (upwelling) and the negative 
values the downward motion (downwelling) . 

From examining the plots of C v and C,. it was 

x y 

determined that when C x and C y vortex lines indicate an upward 
(downward) motion, they were consistent with the positive 
(negative) ^ lines. 

2. Indirect verification of C-Vector vortex lines by SST 

images 

Figure 41 is an AVHRR image of the area for the day 
prior to the cruise (15 May 1991) . It suggests that a plume 
of cold water (8° C) extends offshore of Point Reyes and over 
the continental slope. This area coincides with positive 
values of ^ (upward motions) south-west of Farallones Islands 
(Figure 38). 

South-east of the continental slope and north-west of 
the Pioneer Seamount, there is an area of relatively warm 
water (10° - 11° C) . At this area, has negative values 

(downward motion) . At the south-west side of Pioneer 
Seamount, there are relatively cold waters (9° - 10°) and at 
the same area, ^ has positive values. Therefore, the SST 
images provide useful information concerning possible vertical 
velocity at the ocean surface and are a helpful indirect 
verification of 
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3. Comparison of vertical ADCP velocities with ^ 

function 

Figures 42 and 43 shows the vertical ADCP velocities 
at 20 m and 50 m at the upper part of the figures, and ^ at 
20 m and 50 m at the lower part of these Figures. Notice that 
at 20 m depth and at the south-east of the region, both ^ and 
vertical ADCP velocities are consistent because they show 
negative values. The 'ir has negative values at the center of 
the region, but the vertical ADCP velocities are negative at 
the south-west side of the area. These two regions of 
negative values coincide partially. The rest of the areas 
have positive values. 

At the plots of 50 m, the vertical ADCP velocities 
indicate negative values at the same areas as at 20 m. The ^ 
shows negative values at the east side of the region and at 
the middle and south region. The areas of negative values of 
^ and negative vertical ADCP velocities partially coincide and 
partially differ. 

There is no perfect match. This happens probably 
because the vertical ADCP data are also affected by tides, 
waves or internal motions. Also, the calculation of 'if 
contains some small errors like the round-off error, 
truncation error and machine error. 
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Figure 14. C 1 / f J values at columns 1 and 2 



46 



l (Km 






Figure IS. C x /f J values at columns 3 and 4. 
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Figure 16. C x / f J values at columns 5 and 6 . 
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Figure 17 . C x / f J values at columns 7 and sT 
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Figure 18. C 1 /f J values at columns 9. 
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Figure 19 . C 1 /f J values at columns (from left to right) 9 , 8 
and 7 . 
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Figure 20. C 2 /f J values at sections A and B. 
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Figure 22. C 2 /f J values at section E. 
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Figure 23 . C 2 /f J values at sections (from left to right) E , C 
and A. 
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Figure 25. C 3 /f J values at 50 meters (upper plot] and 200 
meters (lower plot) . 
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Figure 26. C 3 /f' values at three layers (surface, 50 meters 

and 200 meters) . 
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Figure 28. C x /f values ac columns 1 and 2. 
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Figure 29. 
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Figure 34. Cy/ f J values ac sections A and B. 
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Figure 35. C y / f J values at sections C and D. 



67 



CY/F FIT SECTION E 

a 




0.0 10.0 20.0 30.0 40.0 50.0 60.0 



X (KM) 



Figure 36. Cy/ f J values ac section E. 
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Figure 37. C y /f J values at sections (from left 
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Figure 38. t at surface. 
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Figure 40. 'I' at three layers (surface, 50 meters and 200 

meters) . 
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Figure 41. A satellite image of the area (15 May, 1991) . 
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Figure 42. Vertical ADC? velocities at 20 meters (upper plot) 
and i' function at 20 meters (lower plot) . 
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V. SUMMARY, CONCLUSIONS AND RECOMMENDATIONS 



A. SUMMARY 

The purpose of this study was to demonstrate the C-Vector 
method and the potential value of this method in the 
calculation of the ageostrophic circulation using CTD and wind 
data. Sometimes the ageostrophic flow is significant and the 
purpose of this study is to emphasize the importance of 
ageostrophic circulation in waters off California as inferred 
by C-Vector using CTD and wind data. 

Spring is an upwelling season for waters off northern 
California. The SST satellite image of the area shows two 
filaments of cold water. One extends west of the Farallones 
islands and the other stretches offshore of the 1,000 m 
isobath. The cold filaments are associated with positive 
values of the vorticity of C-Vector (equation (29)) . 

The main upwelling extends to a depth of 200 m, and the 
strong C-Vector cells extend from the surface to about 200 m. 
Cells below 200 m are observed to be very weak. This is in 
phase with the Sverdrup assumption that the upwelling water 
rises from moderate depths only, probably less than 200 m 
(Sverdrup and Fleming, 1941) , and the phenomenon represents 
overturning of the upper layer. 
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Some large values of C^, C 2 and C x , C y (horizontal 
components of C-Vector for the ageostrophic and the total 
circulation) are found over the continental slope. This 
probably occurs because of the presence of the continental 
slope and is associated with oceanic processes which occur 
there . 

A shortcoming in the calculations of C-Vector could be the 
uncertainty associated with the determination of a reference 
level. Density differences between stations demonstrate that 
1,400 m is the main representative LNM for the deep waters. 
While the depths where the difference in density between 
stations is zero are not always 1,400 m, there is a good 
indication that at this level currents are very small. 

The Beta- spiral method (Schott and Stommel, 1978) 
overcomes the erroneous "level of no motion" and allows the 
"level" to vary realistically in (x,y). 

B. CONCLUSIONS 

The C-Vector is a good method for calculating the three- 
dimensional pseudo- vorticity fields in the ocean using CTD 
measurements and wind stress data. 

This technique (with the help of 'P) provides some useful 
information of the upward and downward motions in the ocean. 

If the pseudo- vorticity equations can be solved, the 
three-dimensional circulation in the ocean can be obtained. 
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C . RECOMMENDATIONS 

Data from May 1991 gave a picture of the utility of C- 
Vector method and the pseudo-vorticity fields at the Gulf of 
Farallones Islands. 

These data are complicated by such local phenomena as 
upwelling, tides, inertial motions and winds. It is 
recommended that the results of this study should be compared 
with : 

• The pseudo-vorticity fields using data from other sources, 
for example from a current meter array covering the area 
for verification 

• Results of other studies and calculations of C-Vector from 
the same area. 

In the future it is also recommended to use the 13 - spiral 
method in the C-Vector computation. The advantage of the ft - 
spiral method is that it overcomes the problems of the "level 
of no motion" and for this reason it is expected to improve 
the C-Vector method. 

Another way to avoid the uncertainty of the level of no 
motion or the reference level could be the use of satellite 
altimetry data. This data can be used for the determination 
of the geostrophic currents at the ocean surface. Then, an 
integration of the thermal wind equation from surface to a 
certain depth z can provide the geostrophic velocities at 
different depths. 
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